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ABSTRACT 

Recently observed emission lines in the X-ray afterglow of gamma ray bursts suggest 
that iron group elements are either produced in the 7-ray burst or are present nearby. 
If this material is the product of a thermonuclear burn, then such material would be 
expected to be rich in Nickel-56. If the nickel remains partially ionized, this prevents the 
electron capture reaction normally associated with the decay of Nickel-56, dramatically 
increasing the decay timescale. Here we examine the consequences of rapid ejection of a 
fraction of a solar mass of iron group material from the center of a collapsar/hypernova. 
The exact rate of decay then depends on the details of the ionization and therefore the 
ejection process. Future observations of iron, nickel and cobalt lines can be used to 
diagnose the origin of these elements and to better understand the astrophysical site 
of 7-ray bursts . In this model the X-ray lines of these iron-group elements could be 
detected in suspected hypernovae that did not produce an observable 7-ray burst due 
to beaming. 

Subject headings: gamma rays: bursts — supernovae — line: profiles 

1. Introduction 

For the last thirty years 7-ray bursts have been observed, but it has only been in the last 
five years that X-ray, optical, and radio counterparts — the 'afterglows' — have been seen (Costa 
et al. 1997; van Paradijs et al. 1997; Frail et al. 1997). Pinpointing the location of the afterglows 
has allowed some of the bursts to be associated with host galaxies. Absorption lines in the optical 
afterglow, together with redshifts obtained from the host galaxies have allowed redshifts to be 
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obtained for the 7-ray bursts themselves, e.g. (Metzger et al. 1997). Redshifts have been obtained 
between z = 0.35 and z = 4.5, confirming the idea that bursts are cosmological in origin. This 
idea is also supported by their isotropic distribution on the sky, as recorded by over 2700 BATSE 
detections. 

The redshift determinations have led to estimates of isotropic equivalent fluxes in 7 rays, 
ranging from 10 50 ergs to 10 54 ergs. Since the upper bound is approaching the rest mass energy of 
the sun, the astrophysical origin of these sites is fairly restrictive. Less energy is contained in the 
bursts if the 7 rays are strongly beamed, but the beaming can not be stronger than a few degrees, 
based on the lack of detection of 'homeless afterglows' (Meszaros et al. 1999; Dalai et al. 2002; 
Vreeswijk et al. 2002; Vanden Berk et al. 2002). Recent analyses of the break in the afterglow 
signal support the idea that 7-ray bursts are beamed to a few degrees (Frail et al. 2001; Panaitescu 
and Kumar 2001). 

The theoretical problem of what causes the 7-ray bursts can be divided has several parts: 
the fireball model, the energy injection mechanism and the astrophysical site. The mechanism 
for actually producing the 7 rays and the afterglows is described by the relativistic fireball model 
(Rees and Meszaros 1992; Meszaros and Rees 1997; Wijers et al. 1997). After energy is injected into 
material composed mostly of electrons and positrons, it is ejected relativistically. Internal shocks 
in the plasma and external shocks caused by contact with the interstellar medium drive synchotron 
radiation from the electrons in a magnetic field. The spectrum of the afterglow, at least in some 
cases, is fairly well fit to a synchotron spectrum. For recent reviews see van Paradijs et al. 2000; 
Piran 1999. 

This model is fairly independent of the initial energy injection mechanism, as well as the as- 
trophysical site for the bursts. For the injection mechanism, neutrino and antineutrino annihilation 
into electron-positron pairs (e.g., Ruffert 1999), and delivery through a pointing flux such as the 
Blanford-Znajek mechanism, have been discussed (e.g., Brown et al. 2000). 

Two contenders for the site of 7-ray bursts are neutron star - neutron star mergers and 'failed 
supernovae' also known as collapsar or hypernova models. The neutron star models (Eichler et al. 
1989; Mochkovitch et al. 1993) have been shown to be dependent on the way in which general 
relativity is handled in the numerics. Janka (1999), using a post newtonian approximation, finds 
short bursts with at most 10 50 ergs ejected. Salmonson et al. (2001) using a GR approximation that 
is exact in the case of spherical symmetry, find longer timescales and larger energies. The collapsar 
models involve the collapse of a massive star which fails to produce a viable shock for an ordinary 
supernova explosion. Models involve a combination of a rotating black hole and magnetic field 
driving jets along the rotation axis of the black hole (Woosley 1993; Paczyhski 1998; MacFadyen 
and Woosley 1999). Of these two models, the collapsar model has recently gained favor with 
identification, in many cases, of host galaxies associated with the 7-ray burst (Bulik et al. 1999; 
Bloom et al. 1999b). Neutron star binaries would be expected to wander out of galaxies at a rate 
inconsistent with the number of observed associations with host galaxies. Further, in a number of 
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7-ray burst light curves, features strongly resembling a supernova have been found. 

Furthermore, recent observations of K a lines from iron have been identified in the X-ray 
afterglows of four 7-ray bursts (Piro et al. 1999; Antonelli et al. 2000; Piro et al. 2000). This is 
more likely to occur in a supernova like model than in the neutron star merger models. However, 
there are still many unsolved puzzles associated with these lines. The lines are too broad (0.5 
KeV) to be accounted for by thermal broadening, although this may be accounted for by electron 
scattering. Line scattering from reprocesser type models is discussed in, e.g., (Weth et al. 2000; 
McLaughlin et al. 2001; Kallman 2001). If the lines come from material produced by nucleosynthesis 
in a 7-ray burst, then one would expect Nickel-56 as the product of a thermonuclear burn. In most 
cases the lines have been identified as iron, in some cases with independent redshift measurements. 
However, the most recent detection is of nickel, not iron (Reeves et al. 2002). 

Two general classes of models have been proposed for these lines. In one case the iron is already 
present in the system before the 7 rays appear and the lines are produced by the interaction of 
the afterglow with the iron, e.g., the Supranova model (Vietri et al. 2001). In the other case the 
lines are produced by another mechanism and are not directly associated with the afterglow. For 
example, one might imagine that the iron is present in what is left of the star after the burst. An 
accretion disk which has formed around the black hole produces an ionizing flux, that reaches this 
iron, which sits on the surface of a 'funnel ' which has been carved out of the star by the original 
ejection, (cf. Rees and Meszaros 2000). 

Nickel, on the other hand, is likely to be produced by a thermonuclear burn in the silicon 
burning layer of the star. This could be identified by a time dependence in the energy of the 
observed line, due to the decay from nickel to cobalt to iron (McLaughlin et al. 2001). This could 
be principle be detected up to a few days after the event, since the half life for nickel is about six 
days, and the difference in line energy from nickel to cobalt is about half a keV. An alternative 
mechanism suggested for the production of nickel is that first material in the very hot accretion 
disk is completely dissociated into protons and neutrons. A wind from this disk will cool as it 
moves out, and may recombine to an equilibrium nucleus, likely again nickel, (McFadyen 2002). 

In this paper we suggest a variation on the above scenarios. Nickel, produced by a thermonu- 
clear burn in the silicon shell, may be ejected out of the accretion disk surrounding the black hole 
in the collapsar model. This nickel will be partially ionized from below, by the flux leaving the 
accretion funnel. The ionized nickel will decay at a different rate than nonionized nickel, since 
decay of nickel proceeds by electron capture. This will produce a distinctive signature pattern in 
the time dependence of recombination lines. In section 2 we make a plausibility argument to show 
that ejected nickel can remain ionized. We discuss the parameters that determine the optical depth 
of the ejected nickel. We describe the progress of an ionization front through the ejected nickel 
and its dependence on various parameters. In section 3 we discuss the relevant nuclear physics of 
electron capture and beta decay of nickel-56 and cobalt-56. In section 4 we give the fractions of 
these elements as functions of time in this model. In section 5 we give conclusions. 
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2. 



The Model 



In this section we present a model of material ejected from a the center of a 7-ray burst. We 
assume that this occurs by way of interaction between the rotating accretion disk surrounding a 
black hole and the magnetic field. If the 7-ray burst is a collapsar/hypernova, then the material is 
likely to be rich in nickel-56. This will happen if the silicon burning shell with an electron fraction 
of roug hly 0.5 is heated to T > 4 x 10 9 K. This type of burning occurs in ordinary supernova 
explosions. Some of this nucleosynthesis may be deposited in the accretion disk. 

First we discuss the optical depth of the nickel and then we describe the motion of an ionization 
front which travels through the ejected material. 



In this subsection we describe the optical depth of material ejected with fairly high velocity 
{(3 ~ v/c ~ 0.2) from the center of a 7-ray burst. There is an ionizing flux coming from the center 
of the object, presumably from an accretion disk surrounding a black hole. As the material moves 
away from the center it also expands, so we approximate its size as roughly V ~ ilr 2 Ar, where 
r = /3ct, Q, = 27r(l — cos 9) and Ar k, r. Here r is the linear dimension, t is time and c is the speed 
of light. For definiteness, we assume in our calculations that the expansion is homologous, i.e., 
v oc r, and the leading edge moves with v = (3c. This way, the density of the ejecta is independent 
of position, and scales with time as t~ 3 . Since the ejecta are highly supersonic, and the radiation 
force is small except near the inner edge, it is fair to approximate the velocity of each part of the 
ejecta as constant. 

We wish to determine how much of the nickel is fully ionized, as both a function of position and 
time. The first step is to estimate the recombination rate and the ionization rate for the expanding 
material. The ionization rate is given by: 



We use cr(Z) = 2.8 x 10 29 Z 4 /V 3 , for the ionization cross section. The lower limit of integration is 
the binding energy of the final K shell electron, hu = 13. 6Z 2 eV. F{v) is the photon energy flux. 

This flux may have different forms depending on the details of the model. The observed 
ionizing fluxes are much larger than can be supplied by Eddington-limited emission from a disk 
around a black hole of a few solar masses. However, a mechanical and/or poynting wind generated 
by the Blandford-Znajek mechanism or neutrino annihilation can have a much larger energy flux, 
and be converted into heat and radiation when it impacts the funnel wall far from the black hole 
(or from a highly magnetized neutron star), as envisaged by Rees k, Meszaros (2000). In that case, 
the expected spectrum might be more like synchrotron radiation or bremsstrahlung with a very 



2.1. Optical Depth 
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high temperature; in both cases, we may approximate the keV X-ray spectrum with a power law: 
F{v) oc 

This is related to the number flux of the material above the ionization threshold as 

Nj = W r (2) 
Jhu z hv 

We parameterize the time dependence of the ionizing source (in 10 53 photons per second) as 

Ni = Nj^lO 5 ^- 1 (3) 
where time, td is measured in days. After some algebra, we rewrite the ionization rate as 

R ton = 6.5 x (°4) 2 ( 2 4) 2 (> (4) 







where 

P 



f°° En. 1 du 

jvz hv~v* 



r _ J U Z hV V /r\ 
S — POO F.. r ' \°) 



i.e., it is a form factor expressing the ionizing photon number flux with the cross section above the 
edge threshold. For a power-law spectrum F(v) oc we have ( = 

Also necessary is the recombination rate, which we estimate by starting with the recombination 
luminosity (Lang 1980), 

L rec = lQ- 21 n e niT- ll2 VZ A erg s" 1 (6) 

Here temperature is in Kelvin, electron, n e and ion, rii densities in cm -3 , and V is volume. The 
recombination rate per ion can be estimated by dividing by the volume and ion density and also 
by the energy of the photon, hv ss Z 2 13.6eV. Some algebra yields a recombination rate of 



Rrec = ^Fday 1 (7) 



Crec = LI X 10 6 f^Y ( -— ) ( r ) ( "~ if £t ) i ) • (*) 



with time t in days where 

1/2 / M \ /a2\ 3 / 1 - cos 20° \ (Z_\ 3 (A 
V T J VO-lM y V (3 ) V l-cos^ J \28J \56 

Here M is the amount of mass in the ejected material and 9 is the opening angle of the cone. T is 
the temperature of the material. 

With the recombination and ionization rates we can calculate the fraction of nonionized nickel 
in the optically thin region of the material as 
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where the coefficient is 

«--*-(rf=^)(^)(^)(^r (")'(")< w> 

If the luminosity of the ionization source is constant, then the fraction of nonionized nickel decreases 
with time. This may be counterintuitive, since the material gets further from the ionizing source. 
But the ionization rate per atom decreases with flux, i.e., as 1/r 2 oc l/t 2 (since r oc t), whereas 
the recombination rate is proportional to density, i.e., scales as 1/r 3 oc l/t 3 , hence ionization wins. 
However, in case of an emptying accretion disk, or of a spinning-down black hole, the luminosity 
will scale roughly as t -1 , which implies a roughly constant ionized fraction. 

The optical depth of the ionizing photons, r = fnonniazr is an decreasing function of time. 
Therefore an ionization front passes through the material. The details of this ionization front 
determine the ratio of ionized to nonionized nickel, in the material and therefore the rate of nickel 
decay. 



2.2. Ionization Front 

We assume that Nickel 56 is ejected from the accretion disk with fairly high velocity. There is 
a source of ionizing photons, which emits in photons per unit time, 

N = AT /i53 10 53 /^ s" 1 = CsNi/Q day" 1 , (11) 
where TVj is the total number of ions in the material and 

M V 1 (A 



c -= 4xl ° i^J Ur"- (12> 



The equation which describes the passage of the ionization front is: 

N I>53 ^At = nif non R ion nr 3 rlAt + mQ^rj.Arf. (13) 



The photons are either absorbed on their way to the front (second term) or at the front (first term). 
Here r is the linear dimension of the ejected material, and rj = r/r max is the fractional distance 
within that material, while ry. is the position of the ionization front. The solid angle subtended by 
the ejected material is SI = 2n(l — cos 9). We rewrite Eq. 13 as 

3C S _ 3C rec ^ 3 r% drf % 



The solution to this equation is 



C, \ ~ ■■ ■■ 



On 
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as long as (3C rec /2) 1//2 >> during the time the ionization front is passing through the material. 
This is true for all the situations considered here. 

Eq. 15 describes the motion of the ionization front. The time it takes a front to pass through 
the entire mass of nickel depends strongly on a. For example, for the parameters considered here, 
it takes 6 days for the front to pass through the material if a = and 17 days for the front to pass 
through if a = 1. 

We have now described a moving, expanding mass of nickel, which is partially ionized behind 
the front and not ionized ahead of the front. We wish to determine the decay properties of the 
nickel, for that we need the nuclear physics described in the next section. 



3. Electron Capture and Beta Decay 

Here we describe the nuclear physics of the A=56 decay chain. The half life of 56 Ni is 6.075 
days, and it decays almost exclusively by electron capture, 

56 Ni + e~ ^ 56 Co + v e . (16) 

The Q-value of this reaction is 2.136 MeV. The decay proceeds primarily through the 3rd and 4th 
excited states of Cobalt-56 which are at 1.45 Mev and 1.720 Mev with respect to the ground state. 
These states have spin and parity of + and 1 + , while the ground state of nickel is + , so the decay 
proceeds through Gamow- Teller and Fermi transitions. 

If the nickel is completely ionized and in a dilute environment (so free electron capture is 
negligible), then it can only decay by emission of a position, 

56 M ^56 Cq + e + + ^ 

The Q-value for this beta plus decay is 1.114 MeV, which makes the 3rd and 4th excited states in 
cobalt energetically inaccessible. There are three remaining possibilities. The decay can proceed 
through the ground state, 4 + , the first excited state at 0.158 Mev, 3 + or the second excited state 
at 0.970 MeV at 2 + . These are all forbidden and to date no (5 + decay from Nickel-56 has ever been 
seen. The current experimental limits place the half life at greater than 2.9 x 10 years (Sur et al. 
1990). Calculations within the large scale shell model have placed the half life at about 4 x 10 4 years 
(Fisker et al. 1999). For our purposes, we assume that completely ionized nickel is quasi-stable. 
We note that because of the long half life of ionized nickel-56, it has been suggested that it may be 
seen in cosmic ray detectors. 

Cobalt-56 also decays primarily (81%) by electron capture into Iron-56. 

56 Co + e~ ^ 56 Fe + v e . (18) 

Its half life is 77.2 days (Junde 1999). The Q-value, 4.566 MeV, is much higher than in the case 
of nickel, and therefore it has more energetically accessible states in iron into which it may decay. 
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Nineteen percent of the time (Jg+ = 0.19), Cobalt-56 decays by emission of a positron, 

56 Co ^56 Fe + e + + Ve Q Q j 

the majority of which is a 4 + to 4 + transition to the second exited state of iron at 2.085 MeV. This 
decay will still proceed, even if the iron is fully ionized. Therefore, for ionized cobalt, the half life 
will be roughly a factor of 5 higher than for cobalt with two or more electrons. Recent data on 
Cobalt-56 decay can be found in (Junde 1999; Meyer et al. 1990). 



dN Co 


— fnon 


N Ni 


fnon 


dt 


TNi 




dN Fe _ 


= fnon 


N C o 




dt 


TCo 



4. Results: Nickel decay in 7-ray burst ejecta 

We now combine the physics of the previous two sections to describe the motion of the ion- 
ization front with the differing decay times of ionized and nonionized nickel and iron. We show 
that this produces unusual decay patterns that may be observable with an instrument such as the 
Chandra X-ray observatory. 

The rate of change of the relative abundances of Nickel-56, Cobalt-56 and Iron-56 

^ = -fnon^ (20) 

at tni 

^-(l-fnon)I^ (21) 

TCo ' T Co 

+ (l-/n™)#— (22) 

' TCo 

Here TNi = Ti/2 m 2 and T1/2 is the half life of nickel and we use the similar relation for cobalt. 
Before the ionization front arrives, when t < C fa rf {3 ' a \c u = (Crec/Cs) 1 ^ 3 -^ the nonionized 
fraction is essentially unity. After the ionization front passes, the nonionized fraction is given by 
Eq. 9. We approximate the nonionized fractions of nickel, cobalt and iron as the same. In reality, 
there will be slight quantitative differences, since the ionization potential changes by roughly a keV 
from nickel to iron. How this difference translates into f non depends on the shape of the spectrum 
of the ionizing source. 

We can solve the above equations approximately analytically, to get the amount of nickel which 
remains as a function of time and also position in the ionization front. The general expression for 
the fraction of nickel remaining after the front has passed is, 

^(r / ,t)=expl-%^)expf-^ f *—) , t d >C fa rf^ (23) 

TNi I \ T Ni JC fa rf^ C non + t d ) T 

Before the front has passed, 

N(r f ,t) = exp(-t d /T Ni ), t d < Cfsf*. (24) 
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For the special cases of a = and a = 1, Eq. 23 becomes, 

(ft I | \ ~Cnon/T~Ni 

r n0n J r d , t d >C fo r f a = (25) 



^(r/.td) « exp(-C /l r 3/2 /^)exp [-^ (t - C fl r 3 / 2 )] , t d > C fl r 3 / 2 a = 1. (26) 



We show the results of our calculations in Figs. 1 and 2. In Figs la,b we show the fraction 
of nickel as a function of position, rf, in the ejected material, for various times. One can clearly 
see the progress of the ionization front. Ahead of the front there is a constant rate of decay as 
indicated by the horizontal line. Behind the front the material is partially ionized and this slows 
down the rate of decay. In the case where the flux from the accretion disk remains constant, the 
decay is much slower than in the case where the flux decreases as t^ 1 . 

In Fig. 2 we show the total fraction of nickel integrated over position. The lowest curve shows 
the amount of nickel as a function of time in the case of completely nonionized nickel. In the upper 
two curves we show the change with rate of decay of the ionizing source. 

In Fig. 3 we show the amounts of nickel, cobalt and iron as functions of time. We begin with 
relative abundances taken from (Woosley and Weaver 1995), so that the iron group has roughly 
90% Nickel-56 and 10% iron. We then follow the decay at every point in the ejected material as 
the ionization front passes and we integrate over position to look at the total amount of nickel, 
cobalt and iron as functions of time. In these examples, the ionization front passes through the 
ejecta in 1-2 weeks. Strictly speaking that poses a problem for seeing the iron group emission lines 
reported in some afterglows within a day, since it is presumed that we see 7-ray bursts along the 
direction of the funnel and the ejecta. In that case, we cannot see a significant line flux along the 
jet direction until the ionization front has passed through the ejecta. Real life, however, is likely to 
be more complex, with homogeneities in the ejecta and a dependence of properties on the distance 
to the jet axis, allowing the line photons to escape in some directions in some bursts. Alternatively, 
if the ejecta start out hot enough, as may be the case in a neutrino-driven disk wind, they may fully 
never recombine and always be more optically thin. Indeed, the lines are not seen in all afterglows 
and are seen only part of the time in those afterglows where they have been seen. 



5. Conclusions 

The astrophysical origin of 7-ray bursts remains a mystery even three decades after their dis- 
covery. However, technological advances in observing techniques have resulted in a great increase in 
data which is helping narrow down the possibilities. The recent observations of iron recombination 
lines in the X-ray afterglows of 7-ray bursts are an important clue to this problem. 

Here we have discussed the importance of potentially observing the decay timscale of iron 
group lines. The rate of decay will reflect the amount of time newly made Nickel-56 has remained 
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ionized. This information must then be interpreted in the context of the astrophysical site for 7-ray 
bursts . 

We have suggested a model which has Nickel-56 ejected from an accretion disk surrounding a 
black hole. Such a geometry is likely to exist at the center of a 7-ray burst. We have shown that 
the ejected material remains ionized in part due to the ionizing flux coming from the accretion disk. 
The degree of ionization depends on the parameters in the problem such as the ejection velocity 
and importantly, the rate of decay of the flux coming from the accretion disk. We suggest that 
future observations of nickel, iron and cobalt lines several days after the burst are an important test 
of this model. Since the decay of Nickel-56 proceeds by electron capture, observation of a reduced 
rate of nickel decay would be a unique signal of ionization. 

An interesting consequence of our model is that it would allow us to test for the presence of a 
hypernova even if we are observing it off-axis: as the ejected nickel is faster than the normal ejecta, 
it is always outside the general supernova material that explodes more isotropically. This implies 
that the X-ray line flux from the volume behind the ionization front will be scattered more or less 
isotropically, and be observable from all directions. We can estimate the redshift out to which such 
a scattered flux would be visible by noting that the direct flux from the lines is detectable out to 
z ~ 1. As long as the ionization front has not passed through the matter, a fraction of order unity 
of the direct flux will be scattered, and it will be roughly isotropic after scattering. This means 
that the flux towards any off-axis observer will be reduced by a factor 47r/f2j e t ~ 30 relative to 
the line flux for an observer with a viewing direction within the jet. The distance out to which 
we can see the scattered flux is less than the maximum distance for the direct flux by the square 
root of that factor, i.e., corresponds to z ~ 0.2. For much closer hypernova suspects, out to a few 
Mpc distance, even the 7-ray flux from the decaying nuclei could be observed with INTEGRAL. 
In both cases, detection would imply a very unusual event, likely a hypernova, since normally the 
iron-group elements produced in a supernova remain deeply hidden in the ejecta for a long time 
before the overlying layers become optically thin. 

If we were to make observations of delayed nickel decay, it would lend support to the collapsar 
model of 7-ray bursts . Since most of the light from the 7-ray bursts comes in the 7 rays and in 
the afterglow, the central engine is difficult to observe. The type of recombination lines predicted 
here, would, if observed, be a unique window into the center of the 7-ray burst. 

RAMJW is supported in part by NASA (award no. 21098). 

REFERENCES 

Antonelli, L. A., et al. 2000, ApJ 545, L39 

Bloom, J. S., Sigurdsson, S., & Pols, O. R. 1999b, MNRAS 305, 763 



- 11 - 



Brown, G.E., Lee, C.-H., Wijers, R.A.M.J., Lee, H.-K., Israelian, G., Bethe, H.A. 2000, New 
Astron., 5, 191 

Bulik, T., Belczynski, K., k Zbijewski, W. 1999, A&AS 138, 483 

Costa, E., et al. 1997, Nature 387, 783 

Dalai, N., Griest, K., k Pruet, J. 2002, ApJ, 564, 209 

Eichler, D., Livio, M., Piran, T., k Schramm, D. N. 1989, Nature 340, 126 

Fisker, J. Lund, Martinez-Pinedo, G. and Langanke, K. (1999) Eur. J. Phys. 5, 229. 

Frail, D. A., Kulkarni, S. R., Nicastro, L., Feroci, M., k Taylor, G. B. 1997, Nature 389, 261 

Frail, D. A., et al. 2001, ApJ 562, L55 

Janka, H-Th., Eberl, T., Ruffert, M. and Fryer, C, 1999, ApJ 527, L39 
H. Junde, 1999 Nuclear Data Sheets 86, 315 

Kallman, T. R., Meszaros, P., and Rees, M. J., astro-ph/0110654 (2001) 
Krumholz, M. R., Thorsett, S. E., k Harrison, F. A. 1998, ApJ 506, L81 

Lang, K., 'Astrophysical Formulae: A Compendium for the Physicist and Astrophysicist', 1980, 
Springer- Verlog, Berlin. 

MacFadyen, A. I. k Woosley, S. E. 1999, ApJ 524, 262. 

MacFadyen, A. I., talk given at the HEAD 2002, April meeting, Albuquerque, NM. 
McLaughlin, G. C, Wijers, R. A. M. J., Brown, G. E. k Bethe, H. A. 2001, ApJ in press. 
Meszaros, P. k Rees, M. J. 1997, ApJ 476, 232 

Meszaros, P., Rees, M.J., k Wijers, R.A.M.J. 1999, New Astron., 4, 303 

Metzger, M. R., Djorgovski, S. G., Kulkarni, S. R., Steidel, C. C, Adelberger, K. L., Frail, D. A., 
Costa, E., k Frontera, F. 1997, Nature 387, 879 

Meyer, R.A. Fizika(Zagreb) 1990, 22, 153, Wang, G., Warburton, E. K., Alburger, D. A., 
Nucl.Instrum.Methods Phys.Res. A272, 791 (1988), and D.I.Bradley, N.J.Stone, J. Rikovska, 
D. Novakova, J. Ferencei, J. Phys. (London) G12, 115 (1986) 

Mochkovitch, R., Hernanz, M., Isern, J., k Martin, X. 1993, Nature 361, 236 

Paczyhski, B. 1998, ApJ 494, L45 

Panaitescu, A. and Kumar, P., 2001, ApJ 560, L49 



- 12 - 



Piran, T. 1999, Phys. Rep. 314, 575 

Piro, L., et al. 1999, ApJ 514, L73 

Piro, L., et al. 2000, Science 290, 955 

Rees, M. J. k Meszaros, P. 1992, MNRAS 258, L41 

Rees, M. J. k Meszaros, P. 2000, ApJ 545, L73 

Reeves, J. N. et al. 2002, Nature 416, 512 

Ruffert, M., k Janka, H.-T. 1999 A&A, 344, 573 

Salmonson, J. D., Wilson, J. R. k G. J. Mathews 2001, ApJ 553, 471 

Sur, B., Norman, E. B., Lesko, K. T., Browne, E. and Larimer, R. M. 1990 Phys. Rev. C 42, 
573-580 

Van Den Berk et al. astro-ph/0111054 

van Paradijs, J., et al. 1997, Nature 386, 686 

van Paradijs, J., Kouveliotou, C, k Wijers, R. A. M. J. 2000, ARA&A 38, 381 

Vietri, M., Ghisellini, G., Lazzati, D., Fiore, F., k Stella, L. 2001, ApJ 550, L43 

Vreeswijk, P.M., et al. 2002, MNRAS, in preparation 

Weth, C. Meszaros, P., Kallman, T., Rees, M. J., ApJ 534, 581 (2000). 

Wijers, R. A. M. J., Rees, M. J., k Meszaros, P. 1997, MNRAS 288, L51 

Woosley, S. E. 1993, ApJ 405, 273 

Woosley, S. E. k Weaver, T. A. 1995, ApJS 101, 181 



This preprint was prepared with the AAS IATprjX macros v5.0. 



- 13 - 








1 

9 


I I I I 













8 
















o 

"■4— » 





7 






o 








CO 
1— 





6 




\ : 


LL 
















5 
















O 





4 






!z 













3 











2 











1 


i 


i i i 







1 







9 







8 








o 

"■4— » 





7 


o 






CC 





6 


LL 











5 








O 





4 


!z 









3 







2 







1 




0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



Fig. 1. — Fractional abundance of nickel as a function of position of the expanding material. The 
horizontal axis represents relative position, zero nearest the ionizing source and one at the furthest 
edge. The different curves represent different times since the initial ejection. The top curve is after 
one day and the rest are, in decending order 3 days, 5 days, 7 days and 20 days. Ahead of the 
ionization front, the material decays normally, as seen by the horizontal lines. In the left panel the 
number of photons released from the ionizing source is constant as a function of time. In the right 
panel the ionizing source decays inversely with time, iV oc t^ 1 . 



-14- 



o 
o 

03 



CD 

o 



B 
o 




20 30 40 50 
Time (days) 



Fig. 2. — This figure plots the total fractional amount of nickel in the ejected material as a function 
of time. The top curve is for a constant ionizing source while the second is for a source which 
decays inversely with time, N oc tj 1 . The lowest curve shows the fraction of nickel in nonionized 
material. 
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Fig. 3. — This figure plots the total fractional amount of nickel (solid), cobalt (dashed) and iron 
(dotted) as a function of time, for the parameters given in the text and N oc ij 1 . 



